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Abstract 

To probe the universality class of the quantum Hall system at the metal- 
insulator critical point, the multifractality of the wave function ip is studied 
for higher Landau levels, N = 1,2, for various range (a) of random potential. 
We have found that, while the multifractal spectrum f(a) (and consequently 
the fractal dimension) does vary with N, the parabolic form for f(ot) indicative 
of a log-normal distribution of if) persists in higher Landau levels. If we relate 
the multifractality with the scaling of localization via the conformal theory, 
an asymptotic recovery of the single-parameter scaling with increasing a is 
seen, in agreement with Huckestein's irrelevant scaling field argument. 

73.40Hm, 64.60Ak, 71.30+h 



1 



Typeset using REVTgX 



It is an intriguing question to ask whether the transition to the Anderson-localized regime 
in disordered systems is regarded as a phase transition. A suggestion came in 1983 to the 
effect that the wave function at the Anderson transition should be self-similar (fractal) 
on all length scales just as the critical point in a phase transition is characterized by a 
diverging correlation length with the states being scale- invariant .0 To be more precise, the 
self-similarity in critical wave functions extends beyond a single scale transformation, so the 
idea has been subsequently developed into the multifractal analysis.^ 

The criticality of the wave function proposed in general has been analyzed for the quan- 
tum Hall (QH) system in particular,0 where the delocalized states, marginally allowed to 
appear in two dimensions(<i = 2) in the presence of magnetic fields, coalesce into a single 
point (the center of each Landau level) on the energy axis. The multifractal analysis has 
in fact proved to be essential for a full description of the criticality in the QH system as 
recently reviewed by Huckestein.S 

On the other hand, the scaling of localization in QH systems is still some way from 
a full understanding: a finite-size scaling study,! which looks into the scaling of £m/M = 
/(M/£) where £m is the localization length in long strips of width M, indicates that the 
single-parameter scaling holds for the QH system for the Landau index N = 0, while a 
two-parameter scaling seems to be required for higher (N > 1) ones. In an attempt to 
reconcile the apparent lack of universality, Huckesteini has suggested that, if one introduces 
an irrelevant scaling length £ irr to describe £,m/M as a function of two variables, (M/£) and 
(M/£i rr ) ,1 the single-parameter scaling is recovered asymptotically on length scales > £i rr . 
The meaning of the second length scale is not fully clarifiedS except that the existence of two 
parameters is reminiscent of Pruisken's two-parameter scaling that involves the Hall {<J X y) 
and longitudinal (<j xx ) conductivities. 

Here natural questions are: (i) how the multifractality will depend on the Landau index, 
and (ii) can this possibly related to the scaling, which are exactly the purpose of the present 
paper. It has in fact long been recognized&i that the localization length £ wildly grows with 
the Landau index N, where a mapping onto the nonlinear a model gives an estimated £ ~ 



£exp(cjg CBA ) with cxscba oc (N + 1/2) being a xx in the self-consistent Born approximation! 
and i the magnetic length. The N dependence indicates that the scaling has to be discussed 
on much larger length scales for higher levels, which is indeed one motivation for introducing 
the irrelevant length. There, we can control the situation if the irrelevant scaling idea is 
valid: the irrelevant length scale £i rr is shown! to rapidly decrease when the correlation 
length of the random potential a is increased from <C I to ~ £. 

Given this background, it will be intriguing how the multifractality will vary with TV 
and/or o. We can also propose that the multifractal analysis provides a way to probe the 
scaling if we combine the multifractality and the scaling. Janfien§ has actually suggested 
that the scaling amplitude A c = Ywim^^m / M should be related with the multifractal 
spectrum f(a) via A c = [ir(ao — d)}' 1 , where «o is the position of the peak in f(a) against 
a and d = 2 is the spatial dimension of the system. This is obtained by relating the typical 
finite-size scaling variable with the multifractal behavior of the critical wave function through 
the conformal mapping.0'0 

The Hamiltonian of the QH system is given by 

H = Y1 \ NX ) ( N+ -) Hwc(NX\+Y, £ \NX)(NX\V\N'X')(N'X'\ (1) 



2 , 

NX x NX N'X' 

with the Landau wave function \NX), the cyclotron frequency u c , and a random potential 
V. We treat each Landau level separately, assuming that the magnetic field is strong enough 
with \ (NX\V\NX') \ < hu c (N ^ N'). 

Following Andoill! we can express the random potential as an assembly of scatterers 
each having a range d as 

vW = E^-"-" p/2 ^- (2) 

i 

Here the position of the scatterer r, is randomly distributed, which makes the correlation of 
the random potential, (V(r)V(r')) cx (27rcr 2 ) -1 exp(— |r — r'\ 2 /2a 2 ), Gaussian as well. The 
strength, Vi, of the scatterer is assumed to be Vi = ±Vo with an equal probability to make 
the broadened Landau level symmetric. In the short-range limit, the potential becomes an 



assembly of 5-functions, V(r) = J2i Vi^ij — r»). The electronic structure is dominated by 
the dimensionless concentration of scatterers, c = 2n£ 2 ni = rii/riL, where rij is the bare 
concentration and til the degeneracy (per unit area) of each Landau level. 

We consider a QH system of size L under periodic boundary conditions in x and y 
directions. We take system sizes as large as L = (27m^) 1 / 2 £ = 200£ ~ 300£, which are 
considerably larger than those in previous studies.@'0 The Hamiltonian is diagonalized by 
the Lanczos algorithm with no re-orthogonalization imposed in obtaining the eigenvector .0 
This algorithm enables us to obtain 'inner' eigenvalues and eigenvectors accurately with min- 
imal storage requirements for large-scale eigenvalue problems. The computations of desired 
eigenvectors have been carried out in two steps. First, the Hamiltonian is tridiagonalized by 
means of the Lanczos algorithm, and then the tridiagonal matrix is diagonalized in a stan- 
dard way. Second, the eigenvectors are calculated ^from the series of reproduced Lanczos 
vectors and the eigenvectors of the tridiagonal matrix. The size of the tridiagonal matrix is 
taken to be 3Nl ~ 5Nl where Nl is the dimension of the Hamiltonian. This achieves an 
accuracy for normalized eigenvectors ip^ better than || Tiip^ — E^ip^ || = 10~ 10 even at the 
band center, where is the eigenvalue. 

Figure 1 depicts a typical spatial behavior of the wave function |?/;(r)| 2 that has an 
eigenvalue closest to the center of N = 0(a) or iV = 1(b) Landau levels for the delta- 
potentials with the concentration c = 10.0. Sparsely distributed spikes are conspicuous in 
the wave functions when shown on a linear scale, which is the characteristic feature at a 
metal-insulator transition point. This is particularly so for N = 0, while the N > 1 wave 
function is more homogeneously distributed. 

We have then obtained the multifractal spectrum f(a) numerically by the box-counting 
procedure. 00 Namely, we divide the system into N(\) = 1/A 2 boxes where A = £/L is the 
ratio of the box size and the system size. The probability Pj(X) of finding an electron in the 
jth box Bj(X) is 




(3) 
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If we define the normalized g-th moment of the box probability as 



^•(A,g)= N [P / X)]Q , (4) 



we can calculate the singularity strength, 



a{q) = ft hiA ' (5) 

and the multifractal spectrum, 

f(a(q)) = Inn ^ — . (6) 

Equations (6) and (7) yield the multifractal spectrum f(u) in a parameter (q) 
representation. El The lim^o is taken by reducing the box size to the smallest physical size, £, 
with A = £/L. The generalized dimension D(q) is derived as D(q) = [f(oc(q))—qa(q)]/(l — q). 
Specifically, D(2) gives the fractal dimension of the wave function, which has been obtained 
more precisely by the direct determination of the box probabilities.^ 

Figure 2(a) presents the f{a) spectrum of the critical wave functions at the center of 
N — 0, 1, 2 Landau levels, respectively, where a smooth function is obtained in a finite region 
of a. We have checked that the system size is large enough by confirming that /(a)'s for 
two different L = 200£, 300£ agree within statistical errors and that we are indeed looking 
at the critical wave function!! If not, the /(a) spectrum would coalesce into a single point 
(a, f(a)) = (d,d) for extended wave functions, or two points, (a, f(a)) = (0, 0), (oo, d) for 
localized ones in the limit L — » oo. 

A salient feature in Fig. 2 (a) is that (i) reflecting the significant difference seen in Fig.l, 
the f(a) spectra does vary systematically with N, but (ii) for all of N = 0,1,2 f{a) fits 
very well to a parabolic form,El 

ft \ a ( a ~ a o) 2 M\ 

The fact that the critical wave functions at higher Landau levels as well as the lowest one 
possess the parabolic form /(a) indicates a log-normal distribution exp[(lnP— «olnA) 2 /4(ao — 



oQlnA], on the box probability P. Since the position of the peak «o (« for q = 0), and the 
curvature are correlated as seen from Eq. (8), f{a) is characterized by a single exponent 
cto(N), which is a function of N. The value of ao decreases from ao(0) = 2.31 ± 0.02 to 
ao(l) = 2.15±0.02 and further down to a (2) = 2.09±0.01. This is the first key result of this 
paper. We can visualize this in a typical profile of the N = 1 wave function on a logarithmic 
scale in Fig. 3(b). The figure shows that the fluctuation in the wave fucntion, which is 
wild on the logarithmic scale, indeed becomes similar to the N = one (Fig. 3(a)) when a 
linear transformation on logarithmic scale, ln\ip\ — > [(«o(0) — d)/(oto(l) — d)]ln|^>| + const, is 
performed to make the two distributions coincide. The approach of chq to the support of the 
measure (= spatial dimension d) thus means that the amplitude of wave functions becomes 
less wildly distributed. 

Accordingly, the fractal dimension D(2) = 2a(2) — f(a(2)) of the critical wave functions 
approaches to d = 2 as we go from N = 0(D(2) = 1.50 ± 0.06) to iV = 1(D(2) = 1.73 ± 
0.05) and N = 2(D(2) = 1.82 ± 0.05), and, consequently, the density autocorrelation, 
(|^(r)^(r + R) | 2 ) ~ R D ^~ d , becomes longer tailed. We have a smaller value for N = 
compared with that obtained by the lattice modeli or the scaling of the static conductivity]^! 
but the values of D(2) for N > 1 are close to those obtained from the exponents of the 
density-density correlation by Mieck.0 

If we turn to finite-ranged potentials, a > 0, Figure 2(b) shows the f{a) spectrum for 
the critical wave function for N = 1 with o = 0.71£ and c = 3.0. The f(a) spectrum 
continues to be parabolic, but the value of «o significantly increases to «o = 2.32 ±0.03 with 
decreased D(2) = 1.51 ±0.07. The respective values of «o and D(2) almost go back to those 
for N = 0. This is also seen in the profile of the wave function in Fig. 1(c). The resemblance 
to the N = situation provides a sign of the recovery to the A^ = situation. 

We can now look at the scaling amplitude estimated as A c = [7r(a ~ d)]" 1 for N = 1 
as a function of /3 2 = 1 ± {cr/£) 2 in Fig. 4. As the correlation length o is increased, the 
value of A c approaches to a value that is consistent with an argument by Lee et a/.0 on 
the universal conductivity at the transition with A c = l/ln(l ± a/2). This recovery of the 
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single-parameter scaling in the higher Landau level with long-range potential agrees with 
the finite-size scaling result of Huckesteini. 

In conclusion, the multifractal spectrum f(a) at criticality does depend on the Landau 
inded N while a universality still exists in that its parabolic form is preserved, and the 
scaling valid for the lowest Landau level is asymptotically recovered as the potential range 
is increased to ~ £ for N — 1. One will be able to see how the localization exponent v 
behaves via the identification \ jv = «o + a0 by calculating the normalization exponent x 
for higher Landau levels. The N dependence of the fractality will also have an implication 
on the anomalous diffusion coefficient D(q,u) ~ (q/\fuj)~ ri , where rj — d— D{2). 
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Figure captions 



FIG. 1. Squared amplitude of a typical wave function on a linear scale at the center of 
each Landau level for the system of size 300£ for the N = Landau level with a = (top), 
N = 1 with a = (middle), and N — 1 with a = 0.711 (bottom). 

FIG. 2. (a) The /(a) spectra (for \q\ < 4) of the wave functions at the center of the 
N = 0,1,2 Landau levels for short-range random potentials (er = 0). The solid circles 
(squares) represent the system size 200£(300f). 

(b) The same for a finite-ranged random potential, a = 0.7l£((3 2 = 1.5), for iV = 1 
Landau level. 

FIG. 3. Logarithm of the squared amplitude of a typical wave function at the center of 
N = Landau level with a — (top) is compared with the same for iV = 1 with a rescaling 
of the wave function to make its log- normal distribution the same as for N = (bottom). 

FIG. 4. The scaling amplitude A c = [n(a - d)]- 1 for iV = 1 against (3 2 = 1 + (ojl) 2 . The 
horizontal dashed line indicates A c = l/ln(l + \fi). 
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